Critical escape velocity for a charged particle moving around a weakly magnetized 

Schwarzschild black hole 
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We discuss charged particles motion in a spacetime of a weakly magnetized static non-rotating 
black hole. We study under which conditions a charged particle originally revolving around the 
black hole at a circular orbit after being kicked by another particle or photon can escape to infinity. 
We determine the escape velocity for particles at the innermost stable circular orbits and discuss 
the properties of particles moving with near-critical velocity. We show that in a general case such 
a motion is chaotic. 
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I. INTRODUCTION 

The mechanism of jets formation in black holes is one 
of the most intriguing problems of modern astrophysics. 
The matter in the accretion disk and the black hole ro- 
tation might provide sufficient energy to form and sup- 
port the jets. It is plausible that the key role in the 
transfer mechanism of this energy to the jets is played 
by magnetic fields [IHH- Both widely discussed models, 
Blandford-Znajek mechanism 0, H[ aD -d Penrose effect 
for magnetic fields [U, H| , assume that in the vicinity of 
the black hole horizon there exists a sufficiently strong 
magnetic field. This field does not change the black hole 
geometry, but its interaction with charged particles and 
plasma is important. Such black holes are called weakly 
magnetized. Study of such problems in all their complex- 
ity requires 3D numerical simulations of the magnctohy- 
drodynamics (MHD) in a strong gravitational field. (For 
discussion and references, see e.g. Q). 

Quite often when dealing with such a complicated 
problem it might be instructive to consider first its differ- 
ent simplifications, which can be treated either analyti- 
cally, or by integrating ordinary differential equations. 
Motion of a charged particle in a weakly magnetized 
black hole is an important example. In this paper we dis- 
cuss some aspects of this problem. Namely, we assume 
that a non-rotating black hole of mass M is surrounded 
by a static axi-symmetric test magnetic field, which is 
homogeneous at infinity, where it has the strength B. In 
the absence of the magnetic field the geodesic equations 
for the particle arc completely integrable. This property 
is also valid for the motion of the charged particle in a 
weakly magnetized black hole, provided such a particle 
moves in the equatorial plane orthogonal to the direction 
of the magnetic field. This case was studied in details 
in pj-TTcj. A similar problem for weakly charged rotat- 



ing black holes was discussed in [111 ] . The main result of 



this study is that in the presence of the magnetic field 
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the innermost stable circular orbit (ISCO) of a charged 
particle is located close to the black hole horizon. In this 
sense, the action of the magnetic field on a charged parti- 
cle is similar to the action of the black hole rotation on a 
neutral particle. In particular, weakly magnetized black 
holes may play a role of particle accelerators II 2k 
ilar to the fast rotating black holes (see e.g. [141 -1 
the references therein). 

In this paper we consider motion of a charged particle 
in a weakly magnetized black hole out of the equatorial 
plane. We focus on the following problem: Suppose the 
particle revolving around the black hole in the equato- 
rial plane is kicked out of it by another particle or pho- 
ton. Under which conditions can such a particle leave 
the black hole vicinity and escape to infinity? In our 
model the magnetic field far from the black hole is ho- 
mogeneous. Hence, at a far distance from the black hole, 
where its gravitational field is weak and can be neglected, 
the charged particle moves in the homogeneous magnetic 
field. In this case, the corresponding equations are com- 
pletely integrable. However, before the particle reaches 
the spatial infinity, it passes through the region where 
both of the fields, the gravitational and magnetic, affect 
its motion, so that the dynamical system looses its com- 
plete integrability and the motion of the charged particle 
may become chaotic. 

Several examples demonstrating a similar chaotic be- 
haviour are known. For instance, even in the absence 
of the magnetic field chaos arises for the motion of a 
particle when a spherically symmetric metric of a black 
hole is perturbed. For axially symmetric deformation of 
a black hole this was demonstrated in [2(|. There exist 
also several papers discussing a chaotic particle motion in 
the Majumdar-Papapcrtou metric, describing the space- 
time with two or more extremely charged black holes in 
equilibrium [U [22[ ■ In the vicinity of each of these black 
holes the spherically symmetric Rcissncr-Nordstrom met- 
ric is deformed by the common action of the other black 
holes. Numerical analysis of a charged particle trajec- 
tories in the presence of a toroidal magnetic field in the 
Schwarzschild spacetime was analyzed in [23j]. The re- 
sults presented there illustrate that under special condi- 
tions charged particles moving in the vicinity of the black 
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hole under the action of the toroidal magnetic field can 
be ejected to infinity along the axis of symmetry. Even 
in the absence of a black hole, when a charged particle 
moves in a non-uniform magnetic field it often has ei- 
ther trap ped chaotic motion 24j . or chaotic scattering 
behavior [25l [26| . Chaotic motion of a charged particle 
in the Ernst spacetime representing a magnetized black 
hole was analysed in (27j . 

In this paper we demonstrate that the charged parti- 
cle motion near a weakly magnetized black hole is gener- 
ically chaotic. We find the critical escape velocity for 
such a particle required to escape to infinity and dis- 
cuss some properties of the near-critical motion. The 
paper is organized as follows: In Sec. II we discuss our 
model and present an expression for the escape velocity 
for a neutral particle. In Sec. Ill we present the equa- 
tions of motion of a charged particle moving around a 
weakly magnetized Schwarzschild black hole. Section IV 
discusses scattering data for our problem. Dimensionless 
form of the equations and the initial conditions for the 
particle escape from ISCO orbits are given in Sec. V. In 
Sec. VI we give several examples of qualitatively differ- 
ent orbits of a charged particle in the weakly magnetized 
black hole. Basin-boundary analysis of the trajectories is 
applied for an analysis of the charged particle motion in 
Sec. VII. There we demonstrate the chaotic properties 
of the trajectories and determine the fractal dimensions 
in the proper domains. General discussion of the results 
is given in Sec. VIII. In this paper we use the sign con- 
ventions adopted in [28| and units where c = 1 . 



II. ESCAPE VELOCITY FOR A NEUTRAL 
PARTICLE 

Before considering the escape velocity problem for a 
charged particle in a weakly magnetized black hole, let 
us remind the well-known results for a similar problem in 
a simpler case when a particle is neutral and the magnetic 
field is absent. The background Schwarzschild metric is 



ds 2 = -fdt 2 +f- 1 dr 2 
duo 2 = dO 2 + sin 2 6 deb 2 , 



r 2 duj 2 . 



f = l- 



Jl 

r 



(1) 
(2) 



Here r g — 2GM is the gravitational radius of the black 
hole. There exist three commuting integrals of motion. 
Two of them are generated by the Killing vectors 



;(t) 



P 1 3 



di 



>(<t>) 



f 1 3 



(3) 



reflecting invariance with respect to time translations and 
rotations around the symmetry axis. The corresponding 
conserved quantities are the specific energy £ and the 
specific azimuthal angular momentum C z , 



£ 



i r 2 sin 2 I 



(4) 
(5) 



Here m is the mass of the particle, and p M = mu^ are 
its 4-velocity and 4-momentum, respectively. Here and 
in what follows, the overdot denotes the derivative with 
respect to the proper time. The third integral of motion 
is the square of the specific total angular momentum 



2 4A2 



c 2 = 



C 2 



2 2 

r v , 



(6) 



Here we denoted by v ± the following quantity: 



v x = -r0 o . (7) 
Using the normalization condition u 2 = — 1 one obtains 



= £ 2 -U. 



U = f(l+£ 2 /r 2 ). 



(8) 



The motion of the particle is planar. Let this plane 
coincides with the equatorial plane. Then 8 = and the 
effective potential for the radial motion takes the form 



U = U = f(l + £ 2 Jr 2 ) 



(9) 



Consider a particle at the circular orbit r = r , where r 
is the local minimum of the effective potential U. This 
orbit exists for r Q £ (3r s , oo). The corresponding specific 
energy and azimuthal angular momentum are 



£ = 



\/2(r - r g ) 
^r (2r - 3r g ) 



y/2r a - 3r 



(10) 



The ISCO is defined by r a = 3r 9 , which corresponds to 
a convolution point of the effective potential. For ISCO 
we have £isco = 2V2/3 rj 0.943 and |£ z isco| = r g ^/3. 

Suppose now that the particle at a circular orbit col- 
lides with another particle or photon, so that after the 
collision it will move within a new plane tilted with re- 
spect to the original equatorial plane. In a general case, 
all three types of its motion are possible: (i) bounded mo- 
tion, (ii) escape to infinity, and (iii) capture by the black 
hole. The result certainly depends on the details of the 
collision mechanism. For small values of the transferred 
energy and momentum the orbit will be only slightly per- 
turbed. However, for larger values of £ — £ the particle 
can go away from the initial plane and finally can be 
captured by the black hole or escape to infinity. 

In a general case, as a result of the collision, the parti- 
cle will have new integrals of motion: £, C z and C 2 . For 
the case of a neutral particle in the Schwarzschild black 
hole one can easily obtain the conditions of escape in an 
analytical form. To be able to obtain results which allow 
rather simple analysis and presentation we simplify the 
problem and reduce the space of initial data to a one- 
parameter set by imposing the following restrictions: (i) 
the azimuthal angular momentum is not changed, and 
(ii) the initial radial velocity after the collision remains 
the same, f a = 0. Under these restrictions there exists 
only one parameter which determines the motion of the 
particle, namely the new value of its energy. Under these 



3 



conditions, as a result of the collision, the particle ac- 
quires a velocity v ± in the direction orthogonal to the 
equatorial plane [sec (0]. 

After the collision, the total angular momentum and 
the energy of the particle are 



C 2 



2 .2 , r2 



£ = \j£o+ v K r o-r g )/r 



(11) 
(12) 



During the particle motion its polar coordinate 9 changes 
within the interval [arcsin(|£ z |/£),7r — arcsin(|£ z |/£)]. 

As a result of the collision, the total angular momen- 
tum of the particle changes from its original value C 2 to 
C? given by Eq. (fTTj) . The effective potential U defined 
by the new value of C 2 is greater than that before the 
collision and its extrema correspond to r = r± , where 



r± 



£ 2 ±,/£ 2 (£ 2 -3r 2 ) 



and r_ (r+) is a monotonically decreasing (increasing) 
function of C 2 . Thus, immediately after the collision the 
particle is still at the turning point (r = 0), which is 
located between the maximum and minimum of the ef- 
fective potential U. Therefore, the particle will escape to 
infinity if £ > 1, or 



Kl> 



(r» ~ r g ) 



(14) 



In particular, for IS CO we have the escape condition 
> v esc > l/v6; where the last equality corresponds 
to £ = 1. 

The black hole metric (QJ has the evident discrete sym- 
metries 



— s- 2tt - 6 and 9 — >• n - 9 . 



(15) 



These symmetries imply the symmetry of the problem 
with respect to the following transformation: 

C z — > —£ z and v ± — > —v ± . (16) 



III. CHARGED PARTICLES IN A 
MAGNETIZED BLACK HOLE 

We consider now the case of weakly magnetized black 
holes. We assume that a particle has the electric charge 
q and that its motion is affected by the magnetic field in 
the black hole exterior. Namely, we assume that there 
exists a magnetic field in the black hole vicinity which 
is static, axi-symmetric and homogeneous at the spatial 
infinity where it has the strength B. According to the 
procedure given in [Io|, [29| , such a magnetic field can be 
constructed as follows: Because the metric (fTJ) is Ricci 
fiat, the Killing vectors ^ obey the equation 



fr.u 



0. 



(17) 



This equation coincides with the Maxwell equation for a 
4-potential in the Lorenz gauge A^.^ — 0. The special 
choice 



A" 



(18) 



corresponds to a test magnetic field, which is homoge- 
neous at the spatial infinity where it has the strength B. 
The electric 4-potcntial (Tl8|) is invariant with respect to 
the isometries corresponding to the Killing vectors, i.e., 



0. 



(19) 



A magnetic field is defined with respect to an observer 
whose 4-velocity is as follows: 



B^ = -- e ^F Xa u h 



(13) where 



£0123 — +1 



and 



F„ 



A„ 



A,, 



(20) 



det( 5 ^), (21) 



(22) 



For a local observer at rest in the Schwarzschild space- 
time (HJ we have < = f^ 1/2 ^ t) , and Eqs. (JUJ) and @DJ) 
give 



1 



B H = B jl/2 cos q S v _ _ gin 



(23) 



At the spatial infinity the magnetic field is directed along 
the vertical z axis. In what follows, we assume that the 
field is directed upward, therefore, we shall take B > 0. 

In a curved spacetimc a charged particle of mass m 
moving in an external electromagnetic field F^ v obeys 
the equation 



m— = qF%u\ 



(24) 



where D/dr is the covariant derivative defined with re- 
spect to the metric ([TJ , = i M is the particle 4-velocity, 
u^Ufi = — 1, and q is its charge. 



Denote by P M 



qAfj, the generalized 4- 



momcntum of the particle. Then the conserved quan- 
tities corresponding to the symmetries of the problem, 
the specific energy and azimuthal angular momentum, 
are defined as follows: 



= if, 



qB 
2m 



r 2 sin 2 i 



Let us denote 



B 



qB 
2m 



(25) 
(26) 

(27) 
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Using these conserved quantities and the normalization 
of the 4- velocity vector one gets 



i = r i £. 



r 2 sin 2 i 



B. 



(28) 



The p and 9 components of Eq. (|24| give, respectively 

C 2 \ . r g (2C z B-l) 



1 

I 

B 2 



2 :(2r-3r s )(^ + ^ 4 ^ 



2r 2 



(2r 



= - 2 -H 
r 



C 2 cos 9 
r 4 sin 3 9 



-B z sin 9 cos 0. 



(29) 
(30) 



The normalization condition it 2 = — 1 gives the following 
first order equation 



£ 

U eS = f 



r 2 +r 2 f9 2 + U eS , 



l + r 2 sin 2 i 



r 2 sin 2 I 



— B 



(31) 
(32) 



This equation is a constraint. If it is satisfied at the 
initial moment of time, then it is always valid, provided 
the dynamics of r(r) and 9{t) is controlled by Eqs. (f2"9")l 
and (|3"U)) . 

Let us discuss the symmetry properties of Eqs. ()28l) 
(|32p . First of all, these equations arc invariant under the 
transformations 



B4 -S. 



(33) 



Therefore, without the loss of generality, we can consider 
the particle of the positive electric charge (and hence 
B > 0). To consider a particle of a negative charge one 
should apply the transformation (|33|) . In other words, a 
trajectory of a negative charge is related to the positive 
charge trajectory by the transformation C z — > —£ Zl 4> — > 
—(f>. However, after making the choice B > 0, one needs 
to study both cases when C z is positive and negative. 
They are physically different: The change of the sign 
of C z corresponds to the change of the direction of the 
Lorentz force acting on the particle. 

The system (f28|) -([32 j) is also invariant with respect to 
the reflection 9 — > 7r — 9. This transformation preserves 
the initial position of the particle and changes v ± — > —v ± 
[see Eq. ([7])]. Therefore, it is sufficient to consider only 
positive values of the kick velocity v ± . 



IV. FLAT SPACETIME LIMIT AND THE 
SCATTERING DATA 

Suppose a kicked particle escapes to infinity. Let us 
discuss first the asymptotic properties of such escaped 
particle. At the asymptotic infinity the gravitational field 
of the black hole vanishes. Hence, such a charged par- 
ticle is moving in a practically flat spacetimc with ho- 
mogeneous magnetic field. This type of motion is well- 
known (see, e.g., (30|V The corresponding equations can 



be obtained from (|28|) - (|32j) by introducing the cylindrical 
coordinates 



R = r sin 9 , 



r cos 9 , 



(34) 



and taking the limit \z\ — > oo while keeping R finite. The 
asymptotic form of the equations is 

R 2 



t = S, 
z = 0, 



■B, 



R=§-&R, 



1 + 



R 2 + R 2 



Cz_ 

R 2 



B 



(35) 
(36) 
(37) 



A solution to these equations is well-known, it repre- 
sents a helix whose axis is directed along the magnetic 
field B, i.e. parallel to the z axis. If the z component 
of the particle velocity vanishes, its trajectory becomes a 
circle located in a z = const plane. Let R c be the radius 
of the circle, then the particle velocity v is given by 



v = - — LB x R c 

m 



(38) 



Here [a x b] is a vector product of the 3-vectors a and 
b defined in a Euclidean space in the standard way. For 
the uniform magnetic field B defined by (fT8")) we have 



A = -\B x R] 

2 L J 



(39) 



Therefore, the generalized 3-momentum vector of the 
particle P reads 



-q[B x RJ + ±[B x R] , 



(40) 



and the corresponding angular 3-momcntum vector 
about the point of intersection of the z-axis and the 
z = const plane is 



L=[Rx P] 



(41) 



It is directed along the z axis, i.e. L = Ln z , where 
n z is a unit vector which defines positive direction of 
z. If the center of the circle is located on the z axis 
we have R = R c and the z component of the angular 
3-momcntum vector L r reads 



L c = -IbR 2 c <0. 
2 

In our notations (f2~T)) we have 

C C = -R 2 B 1 R c = ^\£ c \/B, 



Thus, for z = the solution to Eqs. (|35J)-(|37]) is 

2B 



R = R c , 



(42) 



(43) 



(44) 



where cf>o is a constant corresponding to t = 0. The 
general solution to Eqs. (|3"5"j) - (f3"T|) can be obtained by 
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coordinate transformation of (|44[) by moving the center 
of the circle to another point R/, = (Rh, 4>h) on the plane 
z = const and boosting the solution in the z direction. 
Accordingly, we derive 



'2 + R? c + 2R h R c cos [fo-^t- <t> h 



= tph + arccos 



R 2 + i? 2 — R 2 
2RhR 



z = zo H ~p. — t . 



(46) 



where zq is a constant corresponding to t = and Vj_oo 
is constant velocity along the z direction. As a result of 
this transformation, the azimuthal angular 3-momcntum 
L z reads 



£ Z = B(R 2 -R 2 C ) = £ C + BR 



h ■ 



(47) 



where Rh is the distance from the z axis to the axis of 
the helix. Thus, one can see that for C z > (C z < 0) 
the z axis is located outside (inside) the circle, while for 
C z = it passes through the circle. 

The energy of the particle does not depend on the lo- 
cation of the circle and can be expressed as follows: 

£ 2 = 1 + i 2 + AB 2 R 2 C = 1 + i 2 + AB 2 R 2 h - 4BC Z . (48) 

If the particle is at rest, we have £ = 1. This corre- 
sponds to z = and R c = 0. According to Eq. (|4"T|) . the 
last equality implies 



£ z = BRi > 0. 



(49) 



Here the equality sign corresponds to Rh = 0, i.e., the 
particle is located on the z axis. If z = and C z < 0, we 
have R c > Rh > 0, and 

£ 2 = 1 + 46 2 i? 2 = 1 + AB 2 R 2 h - ABC Z > 1 . (50) 

Returning back to our main problem of the charged 
particle escape, we can formulate the corresponding scat- 
tering data as the following set: {£oo,R/i}. 



V. DIMENSIONLESS FORM OF THE 
EQUATIONS 

After these remarks we return to our problem. We 
shall integrate numerically the dynamical equations. For 
this purpose we introduce the following dimensionlcss 
quantities: 



<?=-, P=-, i=—, b = Br g . (51) 



The p and 6 components of the dynamical equation 
together with the expression for the energy £ take the 



form 



da 2 



b -(2p-l)sm 2 8, 



(2€6 — 1) £ 2 (2p~3) 



2p 2 



2p 4 sin 2 6 



d 2 9 _ 2dpd6 
la 2 ~ 



cos 8 



p da da p 4 sin 3 i 



£ y 



U eS = 



1 

1 - - 
P 



- b sin 9 cos 9 . 



2 



PiP- 1 ) [fa) +C/ ^' 

' (£ - bp 2 sin 2 6>) 2 



p 2 sin 2 i 



(52) 
(53) 

(54) 

(55) 



The energy of a particle revolving around the black 
hole in a circular orbit of radius p Q at the equatorial plane 
9 = tt/2 is 



i 



i 

Po 



1/2 r 



(±_b P l) 2 
Pi 



1 



1/2 



(56) 



As we did before, in the case of a neutral particle, we 
assume that a kick does not change the particle's az- 
imuthal angular momentum i, but only gives the particle 
the transverse velocity v ± > 0. As a result, the particle 
energy changes from £ D to 



£ = 



£l 



(Po - 1) 
\ 

Po 



1/2 



(57) 



We would like to know whether the particle will escape 
to the asymptotic infinity. To simplify the problem, we 
shall consider a particle initially moving in an ISCO. In 
this case, the parameters I and b are defined by the radius 
of the orbit p Q as follows (see 0] ) : 



±- 



PoiZpo - 1) 1/2 



1/2 



V2 (a p 2 - 9 Po + 3 ± v/(3p -l)(3-p )) 
,_ V2(3- Po )^ 2 

1/2 " 

2 Po (a p 2 - 9 Po + 3 ± ^(3p - 1)(3 - po)) 

(58) 

Here for I > we have p £ (1,3], and for I < we have 
Po £ [(5 + \/l3)/4, 3]. For Po = 3 wc have I = ±V3 and 
b = 0. For such parametrization the magnetic field, b, 
as well as the specific angular momentum I are uniquely 
specified by the radius of ISCO p D , while the only one 
left parameter £ serves to specify the kick. 



VI. TYPES OF THE TRAJECTORIES 

Given the orbit radius p = p Q and the initial energy 
£ > £ of the particle after the kick, we integrate the 
dynamical equations (I52l) - (|53[) numerically. As a result 
of the integration, we can find a trajectory corresponding 
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FIG. 1: Examples of "capture" trajectories. In each case a charged particle is kicked up from its original ISCO at 
p = 2 for I > 0. The energy £ after the kick is (a) 1.12 (b) 1.2 and (c) 1.3. 




FIG. 2: Examples of "capture" trajectories. In each case a charged particle is kicked up from its original ISCO at 
Po = 2.5 for I < 0. The energy £ after the kick is (a) 1.35 (b) 1.475 and (c) 1.525. 



to the given initial conditions. The dynamical equations 
were solved using the built-in Mathematica 8.0 function 
NDSolve. The integral of motion of the system £ [see 
Eqs. (|5"4"|) and ([57]) ] was used to estimate the accuracy 
of the numerical solver. For our calculations the energy 
error was found to be less than 10 -6 . 

Results of the numerical integration show that there 
are three different types of the final particle motion: 

1. The particle is captured by the black hole. 

2. The particle escapes to z — s- — oo. 

3. The particle escapes to z — > +oo. 

The outcome of the motion is considered a capture 
when p reaches 1. It is considered an escape if \z\ 



reaches 10 3 . The maximum computation time was chosen 
a = 10 5 . In escape cases, it was found that the cumula- 
tive error can reach 10~ 2 . The accuracy of the numerical 
solver can be increased to achieve much better accuracy. 
While increasing the accuracy is not a problem when few 
trajectories are to be generated, it can increase the com- 
putation time greatly when the equations of motions are 
intergradcd hundreds of thousands of times as in the case 
of producing basin-boundary plots (see below). However, 
at least in the cases we have studied, increasing the ac- 
curacy of the numerical solver does not change the final 
state of the motion significantly. 

To distinguish the final states of the particle we intro- 
duce an integer number n. It takes the value for the 
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capture case and ±1 for escape to z — > ±00, respectively. 

If one focuses on a single trajectory, one can see that 
there exists a variety of its types with qualitatively dif- 
ferent behavior of the particle within the domain close to 
the black hole. For each of these types of the motion the 
particle may pass a number of times through the equa- 
torial plane. To characterize the dynamical behaviour of 
the particle in this domain we introduce the "winding" 
number w counting how many times the particle crosses 
the equatorial plane before it gets captured or escapes to 
the spatial infinity. The integer number w is, evidently, 
a topological invariant. For our choice of the initial con- 
ditions, the particle after the kick starts its motion from 
the equatorial plane in the positive z— direction. For this 
reason, it is evident, that for n = +1 the winding number 
w is even, while for n = —1 it is odd. In particular, when 
a particle goes to z — > +00 without further crossing the 
equatorial plane, w = 0. 

In rare cases of the computations the particle stays in 
the vicinity of the equatorial plane during the numerical 
lifetime a. It crosses the plane many times and forms a 
compact cloud in the corresponding phase space. How- 
ever, we expect that such a particle eventually falls into 
the black hole or escapes to the infinity. 

Figures HHH illustrate possible "capture" and "escape" 
trajectories of the charged particle near the magnetized 
black hole for both I > at p Q = 2 and I < at p Q = 2.5 
cases. They collect several examples of such trajectories 
for different values of the invariants n and w. 



The fate of a kicked particle was found to be extremely 
sensitive to the initial conditions: even a very tiny change 
of these conditions may drastically modify its global be- 
havior. Such an extreme sensitivity is an indication of 
non-intcgrability of the system and its chaotic nature. 
The dynamical system fl24j) is Hamiltonian and conser- 
vative. Chaos in Hamiltonian systems was extensively 
studied (see, e.g., [3l[ and the references therein). How- 
ever, most of the known tools and results, such as con- 
struction of the Poincare sections and the KAM theorem, 
are related to (quasi) periodic bounded type of motion, 
while in our case the particle's motion can be unbounded. 
For an analysis of our dynamical system we shall use the 
methods developed for the study of chaotic dynamical 
systems which posses several attractors. In our case, we 
have three different types of asymptotic trajectories. For 
such a motion a phenomenon of fractal basin boundaries 
takes place. We shall discuss it in the next section. 



VII. BASIN-BOUNDARY ANALYSIS 

A. Basin-boundary plots 

The chaotic nature of a dynamical system is exhibited 
in its sensitive dependence of trajectories on initial con- 
ditions. A special case is, what is called, the final state 
sensitivity. Such a sensitivity takes place if a dynamical 
system has several coexisting attractors. As a result, for 




given initial conditions, a dynamical trajectory will typ- 
ically converge to one of the attractors. Therefore, there 
must be boundaries of basins of the attractors separating 
the sets of trajectories with different final states. Such 
boundaries are often fractals. The larger is the fractal 
dimension of the basin boundaries the less predictable is 
the behavior of the dynamical system. To measure frac- 
tal dimension one can use the box-counting dimension. 
Details and additional information can be found in, e.g. 
[3ll . l32l ]. This method, in particular, was used in the anal- 
ysis of another dynamical system appearing in General 
Relativity (see, e.g. HU). 

Here we use the basin-boundary method for the analy- 
sis of the asymptotic behavior of the particle trajectories. 
The corresponding plots are presented in Fig. [5] As we 
mentioned, the particle is initially at the ISCO. Plot (a) 
corresponds to the case £ > and plot (b) to the case 
£ < 0. The horizontal axis on these plots shows the 
dimensional parameter p Q of the ISCO, while the verti- 
cal axis shows the dimensionless energy £ of the particle 
after the kick. The white domains in these plots corre- 
spond to the region forbidden for the particle's ISCO. 
The shadowed regions in these plots consist of square 
pixels of the side size p and £ equal to 2.5 x 1CP 3 . We 
used three different colors for these pixels. The color of a 
pixel determines the final outcome of the particle motion. 
These colors are chosen so that dark grey corresponds to 
the particle capture, grey color corresponds to escape to 
z — > +00 (n = 1), while light grey corresponds to escape 



to z — > —00 (n = — 1). 



B. Critical escape energy and velocity 

The general structure of the plots shown in Fig. [5] can 
be described as follows: The dark grey region which ad- 
joints to the white region at the bottom of the plots cor- 
responds to the particle capture. The uniformly grey 
region in the upper-right part of the plots corresponds 
to the particle escape to z — > +00. The winding number 
in this region is w = 0. This region is restricted from 
below by a diffuse domain. The uniform region is sepa- 
rated from the diffuse domain by the line which we call 
the critical escape energy line. Using our numerical re- 
sults we can estimate the critical escape energy of the 
particle describing this line. The approximate analytical 
expression for the case £ > is 

0.115(3.463 -pa) 



f 



(59) 



{p - 1.851) (3.433 -po) 

The maximal relative error of the escape energy is 0.2%. 
Plot[6ja) illustrates the escape energy £ esc . We can ob- 
tain a similar approximate analytical relation for the es- 
cape energy for the £ < case. It reads 

0.4393(3.198 - Po ) 
° sc + (p - 2.105)(3.667 - Po )' [W) 

where the maximal relative error is also 0.2%. This ex- 
pression is illustrated in plot [B^b) . The critical escape 



1-0 1-5 2.0 2.5 3.0 2.2 2.4 2.6 2.8 3.0 

(a) (b) 



FIG. 5: Basin-boundary plots for a charged particle kicked from the ISCO of the given p (horizontal axis) defined 
by the magnetic field b with different energies (vertical axis). The dark grey zones correspond to capture, the grey 
zones correspond to escape toz-> +00, while the light grey zones correspond to escape to z —00. Plot (a): £ > 

and plot (b) £ < 0. The step size for both p and £ is 2.5 x 10~ 3 . 




2.0 2.2 2.4 2.6 2.8 3.0 2.2 2.3 2.4 2.5 2.6 2.7 2.8 2.9 3.0 

(a) (b) 

FIG. 6: Critical escape energy for £ > [plot (a)] and £ < [plot (b)]. The curves are given by the analytical 

expression ([ST?)) for plot (a) and by (|5Dj) for plot (b). 



velocity w^ sc as a function of p Q can be derived using C. Near-critical behavior 

these expressions together with Eqs. (f56|) - (f58|) . Plots 

(a) and (b) in Fig. [7] illustrate the critical escape veloc- For a given po and energies closc but lcss than the crit _ 

ity for £ > and £ < 0, respectively. ical one the final state of the partic i e canno t be strictly 

predicted. The corresponding near-critical domain con- 
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FIG. 7: Critical escape velocity for i > [plot (a)] and for i < [plot (b)]. 



i ■ ii 1 1 1 1 iii ii mil i 

1.90 1.92 1.94 1 96 1.9B 2.00 



(a) A stripe at £ around 1.9 for £ > 0. 




(b) A stripe at £ around 2.5 for £ < 0. 



FIG. 8: Magnified stripes from the fractal region in the 
basin-boundary plots. 



tains the final states of all the three different types. To 
illustrate it more clearly let us consider a small stripe in 
these domains. Magnifications of the stripes are shown 
in Fig. IS] The horizontal axis on these plots shows the 
dimensional parameter p of the ISCO, while the vertical 
axis shows the dimensionless energy £ of the particle af- 
ter the kick which chosen close to 1.9 for I > and close 
to 2.5 for £ < 0. These magnified plots demonstrate a 
linear structure of different regions corresponding to cap- 
ture and both the types of the escape corresponding to 
n = ill. Similar magnified stripes can be constructed for 
different values of the magnification factor. The remark- 
able fact is that each of such plots has similar structure 
which does not depend on the value of the magnifica- 
tion. In other words, the near-critical diffuse domain has 
fractal structure. This fractal structure is a very compli- 
cated Cantor-set like structure, such that a magnification 



of any portion of the fractal region reveals similar pattern 
of the escape and capture regions on a smaller scale and 
it continues ad infinitum. In the fractal regions the wind- 
ing number w corresponding to either escape or capture 
can take different values and generally increases with the 
increasing repetition of the patterns. 



D. Fractal dimension of the near-critical domains 

To get a qualitative measure of the complexity of the 
fractal regions we calculate the box-counting fractal di- 
mension Df, 

, In. Me) , x 

e->0 In - 
e 

where N(e) is the number of squares of the side-length 
e needed to cover a basin boundary. Such squares are 
counted only if they contain at least two different colors. 
The box-counting fractal dimension gives us a quantita- 
tive measure of uncertainty in our numerical computa- 
tions (see, e.g., [3ll HH). Namely, if our current uncer- 
tainty is, say A, and we want to reduce it by a factor 
of 10 _ ™ by improving the precision of our computations, 
then the necessary precision of our computations, say p, 
is given by 

Tl 

log 10 p = log 10 A , (62) 

a 

where 

a = 2-D f , < a < 1 (63) 
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1 2 3 4 5 6 

-lne 



(a) 

FIG. 9: The box-counting dimension. Plot (a): In N(e) vs. 



is the uncertainty exponent. Thus, we need n/a addi- 
tional digits to achieve the desired precision. 

Figure[9]contains plots of In N(e) vs. In (1/ e) for differ- 
ent values of e for the fractal structures shown in Fig. [5] 
The plots illustrate a linear relation for sufficiently small 
e. The fractal dimensions of the two basins-boundaries 
are 

Dwl.60, £>0, (64) 
Dwl.85, £<0. (65) 

The fractal dimension is closer to 2 for the £ < case. 
Thus, the uncertainty exponent a ~ 0.15 is smaller than 
a w 0.40 for the case of £ > 0, and the corresponding 
near-critical domain has more complex fractal structure. 
In this case an increase in precision requires more digits 
in computations. 

E. Additional details of the basin-boundary plots. 

We described the main features of the structure and 
domains in the basin-boundary plots. However, these 
plots contain additional structure which we briefly de- 
scribe now. First of all, let us mention that for t > and 
the values of p in the vicinity of « 1.5 and £ <J 1.2 there 
is an escape lagoon illustrated by the light-grey color. For 
the initial data corresponding to this lagoon the charged 
particle also escapes to infinity but in the direction oppo- 
site to the initial kick, that is with n = — 1. The winding 
number in this region is w = 1. Besides this, there are 
also smaller size light grey regions which correspond to 
different values of the winding number. These regions 




1 2 3 4 5 



-lne 

(b) 

ln(l/e) for £ > 0. Plot (b): In N(e) vs. ln(l/e) for £ < 0. 



form a well- visible set of the light grey stripes located to 
the right of the lagoon. Similar light grey stripes corre- 
sponding to the backscattering to z — ¥ — oo are present 
for £ < but they are much less profound. 



VIII. SUMMARY 

We studied a charged particle motion in the spacctimc 
of a weakly magnetized black hole. We demonstrated 
that the space of its trajectories has rather a rich struc- 
ture. There exist three different types of asymptotic be- 
havior: capture and escape to the asymptotic spatial in- 
finity, z — > ±oo. There also exists a class of escape tra- 
jectories when the charged particle spends a considerable 
time moving in the vicinity of the black hole close to its 
equatorial plane, crossing it again and again. For such 
a particle its escape to infinity has features similar to a 
diffusion process. 

Certainly, our model is rather simplified. We made 
several assumptions that simplify the problem. These 
simplifications are of two different types. First, we chose 
a special type of the magnetic field. In a 'realistic' case 
the magnetic field decreases at far distances. The other 
simplification was the choice of orbits and parameters 
of the 'kicks'. It should be emphasized that the basin- 
boundary plots were constructed for a very special case of 
the kicking mechanism. For other more general types of 
kicks the corresponding basin-boundary plots might be 
modified. However, we expect that the following three 
main features are common: Namely, the motion of the 
kicked particle is mainly chaotic. There is a critical es- 
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cape energy line (surface) and the near-critical domains 
which have fractal structure. It is interesting to confirm 
this by direct numerical calculations. Another interest- 
ing generalization of the problem is an analysis of the 
critical escape phenomenon and structure of near-critical 
domains for rotating black holes in the presence of a mag- 



netic field. 
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